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In this paper, we first study spherically symmetric, stationary vacuum configurations in general 
covariant theory {U{1) extension) of Hoi'ava-Lifshitz gravity with the projectability condition and an 
arbitrary value of the coupling constant A. We obtain all the solutions with the assumed symmetry in 
closed forms. If the gauge field A and the Newtonian prepotential do not directly couple to matter 
fields, the theory is inconsistent with solar system tests for A 7^ 1, no matter how small |A — 1| is. This 
is shown to be true also with the most general ansatz of spherically symmetric (but not necessarily 
stationary) configurations. Therefore, to be consistent with observations, one needs either to find a 
mechanism to restrict A precisely to its relativistic value \gr = 1, or to consider A and/or ip as parts 
of the 4-dimensional metric on which matter fields propagate. In the latter, requiring that the line 
element be invariant not only under the foliation-preserving diffeomorphism but also under the local 
C/(l) transformations, we propose the replacements, N ^ N — v{A — A) /(? and A''* — >■ A^* + NV^(fi, 
where d is a dimensionless coupling constant to be constrained by observations, N and A^* are, 
respectively, the lapse function and shift vector, and A = —if + N^Viip + N{\/iip)'^ /2. With this 
prescription, we show explicitly that the aforementioned solutions are consistent with solar system 
tests for both A = 1 and A 7^ 1, provided that |w — 1| < 10~^. From this result, the physical 
and geometrical interpretations of the fields A and become clear. However, it still remains to be 
understood how to obtain such a prescription from the action principle. 

PACS numbers: 04.60.-m; 98.80.Cq; 98.80.-k; 98.80.Bp 
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I. INTRODUCTION 

Einstein's classical general theory of relativity (GR) is 
consistent with all the experiments and observations car- 
ried out so far fl*]. However, it has been known for a long 
time that GR is not (perturbatively) renormalizable ^ , 
and thus can be considered only as a low energy effec- 
tive theory. Because of the universal coupling of gravity 
to all forms of energy, it is expected that gravity too 
should have a quantum mechanical description. Moti- 
vated by this strong anticipation, quantization of gravi- 
tational fields has been one of the main driving forces in 
theoretical physics in the past decades in a wide range of 
approaches [2]. 

Recently, Hof ava Q proposed a new theory of quantum 
gravity in the framework of quantum field theory. One 
of the essential ingredients of the theory is the inclusion 
of higher-dimensional spatial (but not time) derivative 
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operators, so that the ultraviolet (UV) behavior is dom- 
inated by them and that they render the theory power- 
counting renormalizable. In the infrared (IR) the lower 
dimensional operators take over, presumably providing a 
healthy low energy limit. The exclusion of higher time 
derivative terms prevents ghost instability Q , but breaks 
Lorcntz symmetry, on the other hand. While the break- 
ing of Lorentz symmetry in the matter sector is highly 
restricted by experiments/observations, in the gravita- 
tional sector the restrictions are much weaker [0] (See 
also [7]). The Lorentz breaking and hence the power- 
counting renormalizability are realized by invoking the 
anisotropic scaling between time and space. 



(1.1) 



This is a reminiscent of Lifshitz scalars [8| in condensed 
matter physics, hence the theory is often referred to as 
the Hofava-Lifshitz (HL) gravity. For the theory to be 
power-counting renormalizable, the critical exponent z 
has to be > 3 0, H. Clearly, such a scaling breaks 
explicitly the Lorentz symmetry and thus 4-dimensional 
diffeomorphism invariance. Hofava assumed that it is 
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broken only down to, the invariance under 

t-^t'{t), x^x'{t,x), (1.2) 

the so-called foliation-preserving diffeomorphism, de- 
noted often by Diff(M, T). The basic quantities are 
the lapse function N, the shift vector N'\ and the 3- 
dimensional spatial metric gij, as introduced more than 
50 years ago by Arnowitt, Deser and Misner fio'l, in order 
to quantize gravity. 

Once the general covariance is broken, it immediately 
results in a proliferation of independent coupling con- 
stants [1, [TT|-|l3l| , which could potentially limit the pre- 
dictive power of the theory. To reduce the number of in- 
dependent coupling constants, Hofava introduced two in- 
dependent conditions, the projectability and the detailed 
balance The former requires that the lapse function 

be a function of t only, 



N - Nit), 



(1.3) 



while the latter requires that the gravitational potential 
should be obtained from a superpotential Wg , where Wg 
is given by an integral of the gravitational Chern-Simons 
term over a 3-dimensional space, Wg ^ J^ujsiT). With 
these two conditions, the general action contains only five 
independent coupling constants. The detailed balance 
condition has several remarkable features [l[ll[l|. For 
example, it is in the same spirit of the AdS / CFT corre- 
spondence [l^, where a string theory including gravity 
defined on one space is equivalent to a quantum field the- 
ory without gravity defined on the conformal boundary 
of this space, which has one or more lower dimension(s). 
Yet, in the non-equilibrium thermodynamics, the coun- 
terpart of the superpotential Wg plays the role of en- 
tropy, while 6Wg/Sgij represents the corresponding en- 
tropic force [l^. This might shed further lights on the 
nature of gravitational forces, as proposed recently by 
Verlinde [13|- For details, we refer readers to Hofava's 
original paper [3] , as well as his review article [l^ . 

When applying the theory to cosmology, various re- 
markable features were found. (See for a review.) In 
particular, the higher-order spatial curvature terms can 
give rise to a bouncing universe [13] , may ameliorate the 
flatness problem 21 1 and lead to caustic avoidance [23 |: 
the anisotropic scaling provides a solution to the horizon 
problem and generation of scale-invariant perturbations 
without inflation (23| , a new mechanism for generation of 
primordial magnetic seed field [IJl, and also a modifica- 
tion of the spectrum of gravitational wave background via 
a peculiar scaling of radiation energy density (25| : with 
the projectability condition, the lack of a local Hamilto- 
nian constraint leads to "dark matter as an integration 
constant" the dark sector can also have its purely 
geometric origins [l^l; in the parity- violating version of 
the theory, circularly polarized gravitational waves can 
also be generated in the early universe [l^; and so on. 

Despite of all the above remarkable features, it was 
found that the projectability condition leads to several 



undesirable properties, including infrared instability d, 
^ and strong coupling [s^, HH ^ . All these properties 
are closely related to the existence of a spin-0 graviton 

[Sill. 

It should be noted, however, that the infrared insta- 
bility does not show up under a certain condition (Toj 
and that the strong coupling is not necessarily a prob- 
lem if nonlinear effects help recovering GR at low energy. 
Of course, the strong coupling implies that the naive 
perturbative expansion breaks down and that a proper 
non-perturbative treatment is needed. In general, non- 
perturbative analysis is not easy to perform in practice. 
Nonetheless, in some simplified situations, fully nonlin- 
ear analysises were already performed, showing that the 
A — ?► 1 limit of the theory is continuous and that GR is 
recovered in a non-perturbative fashion. Such examples 
include spherically symmetric, stationary, vacuum con- 
figurations [19|, a class of exact cosmological solutions 
31 [ and nonlinear superhorizon perturbations .35., ^3Q]. 
The non-perturbative recovery of GR, explicitly shown 
in those examples, may be considered as an analogue of 
the Vainshtein effect |37| . 

Although the existence of spin-0 graviton after all may 
not be a problem due to the analogue of the Vainshtein 
effect, it is interesting and certainly important to seek 
another possible way out. Motivated by this, Hofava 
and Melby- Thompson (HMT) [H] extended the symme- 
try ([1^2)) to include a local C/(l), 



U{1) K Diff(M, T). 



(1.4) 



With this enlarged symmetry, the spin-0 graviton is elim- 
inated [1^ , and the theory has the same number of 
propagating degrees of freedom as GR. This was initially 
done in the special case with A = 1, and was soon gen- 
eralized to the case with any A Even with A ^ 1, 
the spin-0 graviton is still eliminated [10, |4l| . When ap- 
plying it to cosmology, various interesting results were 
found [4^ . In particular, the Friedmann-Robterson- 
Walker (FRW) universe is necessarily fiat in such a setup, 
provided that the coupling of the U(l) field to a scalar 
matter field is described by the recipe given in [ioj . 

In this paper, we shall consider two important issues in 
the general covariant theory of the HL gravity with the 
projectability condition (|1.3p and an arbitrary coupling 
constant A [3l,|4^|4l|: (i) the solar system tests; and (ii) 



^ Note that even without the_projectabiHty condition the theory is 
still strongly coupled 12, ,331, although instability can be avoided 
by inclusion of the term a^a' [Till , where ai = N^i/N. On the 
other hand, as mentioned above, abandoning the projectabil- 
ity condition results in a proliferation of independent coupling 
constants. To render this problem, Zhu, Wu, Wang, and Shu 
recently introduced a local U{1) symmetry (See Ea.l ll.4l l given 
below), in addition to the detailed balance condition ,3S] . In 
order to have a healthy IR limit, however, they found that the 
latter has to be broken softly by adding all the low dimensional 
relevant terms. Even with these terms, the number of the inde- 
pendently coupling constants is reduced to 15. 
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the physical and geometrical interpretations of the gauge 
field A and Newtonian pre-potential ip. Specifically, after 
giving a brief introduction to the theory in Sec. II, we 
present all the spherically symmetric, stationary, vacuum 
solutions of the theory in closed forms in Sec. III. In Sec. 
IV, we consider the solar system tests by not taking A 
and if as parts of the low energy 4-dimensional metric on 
which matter fields propagate, and find that theory is not 
consistent with observations as long as A is not precisely 
equal to one, however small |A — 1| is. In Sec. V, we 
further study the limit A — >■ 1 without assuming that the 
configuration is stationary. We find that, although the 
limit exists, it does not reduces to the Schwarzschild ge- 
ometry. These resuts in Sec. IV and V strongly suggest 
that, in order for the theory to be consistent with the so- 
lar system tests, A and/or ip should enter the low-energy 
4-dimensional metric. Thus, in Sec. VI, by requiring that 
the line element ds^ be gauge-invariant not only under 
Diff(A/, J^), but also under the U{1) transformations, we 
propose that it should take the form. 



ds^ = -Af^c^dt^ + 5y {dx' +N'dt) {dx^ +N^dt) 



where 



(1.5) 



N = N-—{A-A), A/"* = TV'' + iVV V, 

A = -ip + N'V,^+]^N{V,ipf, (1.6) 

where w is a dimensionless coupling constant to be con- 
strained by experiments/observations, and subjected to 
radiative corrections. denotes the covariant derivative 
with respect to the 3-metric gij . With such replacements, 
in this section we show explicitly that the resulted met- 
rics are consistent with observations for both A = 1 and 
A ^ 1. With these replacements, one also sees clearly the 
physical and geometric meanings of A and Lp. Our main 
results are summarized and discussed in Sec. VII. 

Note that solar system tests were studied in other ver- 
sions of the HL theory previously [43l |. However, to our 
best knowledge, in the current paper it is the first time to 
consider the problem in the general covariant theory of 
the HL gravity with the projectability condition and an 
arbitrary coupling constant A, while the case with A = 1 
was studied in |44l |. 

In addition, all the high-order derivative terms of cur- 
vatures are negligible in the IR. Then, test particles move 
along geodesies, as shown explicitly in El] by using 
optical geometric approximations. Therefore, to have 
a consistent treatment, when we consider solar system 
tests, we ignore all the corrections from these high-order 
terms. 



II. GENERAL COVARIANT HL THEORY 

To realize the enlarged symmetry (|1.4I) , HMT observed 
that the linearized (minimal) HL theory has a global U(l) 



symmetry for A = 1. This symmetry can be promoted 
to a local one by introducing a gauge field A, with which 
it was found that the scalar degree of freedom is elimi- 
nated [38| . When they lifted it to a full nonlinear theory, 
HMT found that the realization of the symmetry (|1.4p 
requires introduction of an auxiliary scalar field (p, which 
was referred to as the "Newtonian prepotential." Under 
the local U(l), both A and ip transform as, 



5aA = a - N'-^/ia, Sa^p = 



(2.1) 



while the lapse function N, the shift vector iV* and the 
3-metric gij transform as, 



0, 



(2.2) 



where a denotes the U(l) generator, and a = da/dt. 

Under the coordinate transformations (II. 2p . ip trans- 
forms as a scalar, while A transforms as a vector under 
the time reparametrizations t — !■ f{t'), and as a scalar 
under the spatial transformations x — > C(t',x'), namely, 



6A - Cd,A + fA + fA, 
6ip = fip + Cdi(p. 



(2.3) 



The metric components, N, and gij, on the other 
hand, transform as 



SN = C^kN + Nf + Nf, 



5 gij = VjCj + VjCi + fg 



(2.4) 



under ([TT 

The HMT model was initially constructed in the case 
A = 1, and it was soon found that it can be generalized 
to the case with an arbitrary A [40|, in which the spin-0 
gravitons are also eliminated [40l.l4ll| , so the gravitational 
sector has the same degree of freedom as that in GR, i.e., 
only massless spin- 2 gravitons exist. 

For any given coupling constant A, the total action can 
be written as fH-Iill, 



5* 



dtd-'xNy^(^CK - Cv 

where g = det gij , and 
Ck 



(2.5) 



r — 



2ii:.y + v,Vj(p 



'pg''[ 

-f2A, 
(l- A)[(vV)V2ifVV 



R 



(2.' 



Here is a coupling constant, and the Ricci and Rie- 
mann terms all refer to the three- metric gij , and 

1 



Kr. 



2N 



1 



-gijR + Ag^y . 



(2.7) 
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Cm is the matter Lagrangian density, which in gen- 
eral is a function of ah the dynamical variables, C/(l) 
gauge field, and the Newtonian prepotential, i.e., Cm = 
Cm{N, Ni, gij, ip, A; x)j where x denotes collectively 
the matter fields. Cy is an arbitrary DifF(I])-invariant 
local scalar functional built out of the spatial metric, its 
Riemann tensor and spatial covariant derivatives, with- 
out the use of time derivatives. 

Note the difference between the notations used here 
and the ones used in [H, 113] ^. In this paper, without 
further explanations, we shall use directly the notations 
and conventions defined in [47'| and {s^. 

In [ill , by assuming that the highest order derivatives 
are six, the minimum in order to have the theory to be 
power-counting renormalizable [3,'9'|, and that the theory 
preserves the parity, the most general form of Cy was 
constructed and is given by. 



The corresponding timelike Killing vector is = dt- In 
the diagonal case, we have fj, — —oo. With the gauge 
freedom of the local U(l) symmetry, without loss of the 
generality, we can always fix the gauge by setting 



ip = 0. 



(3.3) 



Cv = 



1 



C^go + 9iR+^ (ff2-R' + g3R^JR'') 
+1 (g4i?' + g^R R^,R'' + geR^RiR'l 



Then, 


we find that 










= 0, C^ = C\ - 


= 0, (3.4) 


cind 
















2;/ 

= — -^r-Jj^ + e-^'' 


Tv' - (1 - e^-)' 






= [2A + (A- 

r '■ 

+re^'-•' [2A - 1 


+ Ar/i'] 




Ck 


= -^[4^- 


2 + 4ArAi' + r2(A- l)(Ai')'] , 



Ca = 



2A 



(1 - irv') + (Agr^ - l) 



(3.5) 



ff7(Vi?)2+58(V,i?,fc)(V'i?J'^) , (2.8) where =(5W+sin2 



where the coupling constants (7s(s = 0,l,2,...8) are all 
dimensionless, and 



= uj I and A — A(r). The expres- 

sion for Cy is very complicated and shall not be given 
explicitly here. 

In the vacuum case, we have 



(2.9) 



— Ja — Jip — 0, Ji — 0, Ti- 



0. 



is the cosmological constant. The relativistic limit in the 
IR requires 



'ip — o/j — U, /.y U. (3.6) 

Then, the Hamiltonian and momentum constraints (jA.ip 
and (|A.2[) reduce, respectively, to 



51 = -1, C 



1 



(2.10) 



j r^e''{CK + Cy)dr = Q, 
a{r)h" + h{r)h' + c{r)h^Q, 



(3.7) 
(3.8) 



Then, the corresponding field equations are given in 
Appendix A. 



III. SPHERICAL VACUUM SOLUTIONS 

Spherically symmetric static vacuum spacetimes with 
projectability condition in the HMT setup were studied 
systematically in ;44i446l lioj . In particular, the ADM 
quantities can be cast in the form [5Cl ISl] , 



where 

a{r) = 

b{r) = 

c(r) = 
with 



(l-A)r2/2, 
\{l-\)rf {4f-r,n, 



(1-A) p (//"-/'')+4r 

h{r) = e^'-r 



rff, (3.9) 
(3.10) 



TV = 1, N'd, = e^'-''^r, 
g.jdx'dx^ = e^^dr^ + r^d'^n, 



(3.1) 



in the spherical coordinates = {r,9,(t)), where d^il 
de'^ + siv? edd)^. and 



(3.2) 



In particular, we have Ki 



j^HMT g_ _ qHMT^ where quantities with the super-indice 



"HMT" are those used in |38l |4l 



f{r) = e-' 
Equation (|A.4p reads, 

(1 - \){rh"' 4/i") + d{r)h' + e{r)h = 0, (3.11) 

where 

dir) = i;?p{4/' + 3(l-A)r2(/')' 

+ 4/[A,r2-l-(l-A)r2/"]}, 

e{r) - ^{3(A-l)r2(/')^ + //'[2(l-rX) 

+ (l-A)r(4/' + 5r/")] +2/2[(2A-l)/' 
+ 4rA,-(l-A)r(2/"-fr/"')]}. (3.12) 
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Equation (jA.5[) . on the other hand, yields, 

(r/)' - (1 - A,r2) = 0, 
while the dynamical equations (|A.7p read 
A \' Gir) 



/1/2 J 2r/3/2 



(3.13) 



(3.14) 



2r^fA" + r (2/ + rf ) A' + r (/' + 2A3r) A + H{r) = 0, 

(3.15) 

where G(r) and H{r) are defined in Eq. ()B.ip . 

It should be noted that not all of the above equa- 
tions are independent. In fact, Eq. (|3.1ip can be obtained 
from Eqs.dMj and (0311), while Eq. (|3TT5l) can be ob- 
tained from Eqs. dXTil) . and ([XT^ . Therefore, in 
the present case there are only three independent differen- 
tial equations, li3.^] . h3.13\) . and Eqs. \3.H^ , for the three 
unknowns, (/, h, A)^ . In particular, from Eg. (13. 131) . we 
find that the general solution for / is given by 



f{r) = 1 - — - 3 Aff^ 



(3.16) 



where B is an integration constant. 

Note that the momentum constraint p.Sp is a linear 
second-order ordinary differential equation for h{r), and 
in principle one can integrate it to find h{r) for the gen- 
eral solution /(r) given above. Once h{r) is found, one 
can integrate Eq. (|3.14p to obtain A, 



A(r)=/V2(^)(Ao-^ 



G{r)dr 
rfy^{r) 



(3.17) 



where Aq is an integration constant, and G{r) is given 
by Eq.lEH). 

On the other hand, for the general solution (I3.16p . the 
potential Cv defined by Eq. p.Sp is given by 45], 



Cv = 2A 
where 



36r9C'* 



{ao + air + a2r^ + asr^) , (3.18) 



-216B3 (56 
ai = 32405^38, 



3058) 



a2 



= 216^2 



(255 + 256 - 558) Ag + 53 



as = 8A, 



4 (954 + 355 + 56) + 6C' (352 + 53) Ag 



(3.19) 



All the solutions with A = 1 were found in j4J, |49[ , and 
their global structures were systematically studied in [4^ . 



^ Certainly, such obtained solutions must satisfy the global con- 
straint 1 13.71 1. 



Therefore, in the rest of this section, we consider only the 
case where A ^ 1. The case A ^ 1 was also studied in 
[5^ . but only approximate solutions were found. 

When A 7^ 1, a particular solution of Eg. p. 81) is h{r) = 
0. Then, from Egs. (|B.l|) and ([3TT6)) we can see that A{r) 
now is independent of A, and the corresponding solutions 
will be the same as those given in the case X — l,h — 
0, / 7^ [i^,!!^. Therefore, in the following, we consider 
only the case where /i(r) 7^ 0. It is found convenient to 
consider the four cases, Ag — = B; Ag = 0, B ^ 
0; Ag 7^ 0, -8 = 0; and BAg ^ 0, separately. 



1. Ag=0 = B 



In this case, the momentum constraint p.8p reduces to 
r^h" + 2rh' ~2h^0, (3.20) 
which has the general solution. 



h{r)=Cir + ^, 



(3.21) 



where Ci and C2 are two integration constants. Inserting 
it into Eg. (PT7|) . we find that 

A{r) =A,-^ + - [3(1 - 3X)C! + 2A] r'. (3.22) 



2. A9 = 0, B / 

When Ag = and B 0, the momentum constraint 
p.8|) reduces to 

r (r - 2B) h" + (2r - hB) h' 
2 fr^ -5Br + 5B'^ 



r \ r-2B 



Bw]h^O, (3.23) 



where -cu = 1/(A — 1). Note that when A = 1, we must 
have B = 0, and Eg. (|3.23l) is identically satisfied for any 
h{r), as noticed previously. When A 7^ 1, setting x — 
r/{2B), h{x) = ho{r)hi{x), Eq.(|32Sl) takes the form. 



j2 L JL 

x{l-x)^+p{x)-l~q{x)hi=0, (3.24) 
ax^ ax 



where 



Pi-) - 2x(l-.)^ + ^i^, 

Hq Ax — 5 h'n 
q{x) = x{x - 1)-^ + 



ho 2 ho 
\x^ - lOx + 5 w 



(3.25) 



2x(l — a;) X 
Assuming that Eg. p.24p takes the form of the hyperge- 



ometric differential equation, 

dx"^ 



x(\ - -I- [c - (a + 6 + 1)2:] — - ahhi = 0, 



(3.26) 
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where a, b and c are constants, from Ea. p.25p we find 

2x{l~x)-^ + -—-=c-ia + b+l)x, (3.27) 
ho 2 

h'J, 4x-5h'r, - lOx + 5 



2a;(l - x) 



— ~ ~ab. 

X 

(3.28) 

Ea. p.27p has the solution, 

ho{x) - {X - l)K3+2a+26-2c)^i(2c-5)^ (3 29) 

for which Eq. p.28p is satisfied identically, provided that 

(a -6)2 -9 = 0, (3.30) 
2ab-c{a + b+l) + 2{5 + w) = 0, (3.31) 
- 8c - 45 - mm = 0. (3.32) 

A solution of the above equations is given by 

a = J(7 + Ao), &=i(Ao-5), 
1 



c = -(2 + Ao), Ao = V49 + 16m. 
Thus, the general solution of h{r) takes the form, 



(3.33) 



Hr) = hoir) 



ai F (a, b; c; a;) 



7 — An —5 — An 2 — An 



+a2 x'^^°F 



(3.34) 

where oi and a2 are two integration (possibly complex) 
constants, and F{a, b; c; z) is the hypergeometric func- 
tion dl] with F{a,b;c;0) = 1. Inserting Eq. (P35)) into 
Eq.(I32Hl), we find that 



ho{r) = 



\2B J \2BJ 
On the other hand, for Ag = Eq. (|3.17p becomes 



(3.35) 



A{r) 



1 Mo 

r 



P{r) 



1 _ 2S 



-.dr , (3.36) 



where 
P{r) 



[(9-25A)B2 + 2(l-2A)r 



4(2B- r)3 
+4Br(5A - 2)] ~ 2r{2B - r) [(5A - 1)B 

-2r\]hh' - {r-2Bf 

2B^ 



AB 12B^ , 
r t' 



-2556 - 2058) 



2^2 



-53 



2Ar2-r2(l-A)(/i') 



(7255 + 8156 - 6358) 



/\2 



(3.37) 



The Hamiltonian constraint (|3.7p now reads, 
r2 (£v+£K)^ir 



1 _ 2B 



0, 



(3.38) 



where Cy is given by Eg. (13. 181) with A^ = 0, and 
(2 B^ \ , . 2Bhh 



r2 (r - 2S)2r2 ) 2Br - r"^ 



5B - 2r 
r{2B - r)' 



+{h'y -\[h' + 



h] . (3.39) 



3. 5 = 0, Ag / 

When B — 0, Ag ^ 0, the momentum constraint (|3.8p 
reduces to 



h" + -[l 

r 



kgr"^ - 3 



h' 



18(A - 1) - 3(5A - 7)Agr2 (A - 3)A2r 



(A-l)r2(Agr2-3)2 



' h = 0. 



(3.40) 



Note that when A = 1 we have Agh = 0. For A 7^ 1, 
Eg. (13.401) has the general solution. 



where 61 and &2 are constants, and 



z = g Agr , Ai = 
Then, Eg. dXTTl) becomes, 

A(r) = ^3 - A<,r2 |^Ao - y 



A-3 
A- 1' 



(3.41) 



(3.42) 



=dr , (3.43) 



where 
Qir) 



12(3 - Agr2)3r 1 
+(A<,r2 - 4)A<,r2 



A(Agr2 - 6)2 - 3 



/i2 + 18(3-Agr2)7 



6A 



(A + l)Agr2 /i/i' + (3 - Agr2)2r2 



6A„ 



4A| 



-I8A + 9(A - l){h'f + ^(352 + 53) 

g^3 

+^(954 + 355 +56) 



(3.44) 
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The Hamiltonian constraint p.7|) . on the other hand, 
takes the form, 



[Cv + Ck) dr 



= 0, 



(3.45) 



where Ck takes the same form of Eq. (13.391) but now with 
h{r) given by Eg. (13.411) and Cv given by. 



Cv 



"3 



36C'*' 

as can be seen from Eqs. (|3.18|) and (|3.19p . 

I B / 0, Aj / 

In this case, Eqs. p.Sp becomes, 
„3^ uii 



(3.46) 



(6B - 3r + kgr^) h" + (15B - 6r + A^r^) /i' 
/6B-2Agr3 18(5^2 -5Sr + r2) 
~ r(A-l) ^ r(6B - 3r + A^r^) 



3Ag(16B-5r)+A2r6^ 
r(6B - 3r + Agr3) 



/i = 0. (3.47) 



Setting 



r = 2Bx, Ao = -B^Ag, (3.48) 



and /i(r) = h(){x)hi{x), where 



/■ 2c - 5 - 2a; [a + - 1 + (Ao - e)x^] 

'^^("^ = ""Py 4x(l-x + Ao) 

(3.49) 

we find that Eq. p.47p reduces to 

x{l~ x + Kox^)h'^ +[c- {a + b+l)x + ex^] h'^ 



{ah + kx^)ho^O, (3.50) 



but now with 



Ao + 7 , Ao - 5 

Ao + 2 Ao 
c = — ^ — , e== — (Ao + 5), 

OA 

k = -^(Ao + 7 + 8n7). (3.51) 
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When c is not an integral, expanding ho{x) in the form, 

oo oo 

/io = Ai^a,a;^+A2^M'+'^", (3.52) 

1=1 4=1 

where Ai and are two constants, we find that in terms 
of the two arbitrary constants aq and bo, the coefficients 



ai and hi {i ^ 0) are given by, 
ai 



ah 

— Oo, 

c 



a2 = 

bi - 
b2 -- 

and 



ab{a + b + ab + 1) 
2c(c4- 1) 
A: 



-ao, 



3(c + 2) 

ab{a + b + ab + l)(2a + 2fe + afe + 4) 



ao, 



6c(c+l)(c + 2) 
(A: - e)ab {a + b + ab + l){2a + 2b + ab + 4) 



4c(c + 3) 
x{3a + 3b + ab + 9)ab + 



24c(c+l)(c + 2)(c + 3) 
k{3a + 3b + ab + 9) 
12(c+2)(c + 3) 



ao, 



afe+ (1 - c)(a + fe - c + 1) 
2-c 

afe+ (1 - c)(a + 5 - c + 1) 

2(2-c)(3-c) 
X [a6 + (2 - c)(a + & - c + 2)] 5o, 



(3.53) 



a6+(j-l)(a + 6 + j-l) 

0-7 — TTT : ^ J — 1 



j(i-l + c) 
(j-3)(jAo-4Ao + e)-fc 
j(j - 1 + c) 



aj-3, i > 5 



(j-2-c)[(j-3-c)Ao + e]-fc ^ 

= — — ^ 0j_3 



jij-3 + c) 
afc + ( j — c) (a + 6 + j — c) 
j(j-3 + c) 



&J-1, j > 3 

(3.54) 



Note that one can always set Ai = A2 = l,hy redefining 
the two arbitrary constants oo and &o- 

On the other hand, when c = 1 + m (where m is 
an integral), we can use the Frobcnius method to solve 
Eg. (13. 501) . Let us first write hQ{x) in the form. 



ho = Ai ^ aiX^ + A 



2X 



i=0 



In X aiX^ + bi 



1=0 



(3.55) 



where 



c — (a + 6 + l)a; + ex^ 
1- x + Aqx^ 

—abx — kx^ 



1 - a; + Aox^ 



2=0 

00 

^^d,x\ 



(3.56) 



then, we can obtain the coefficients a^, at and bi in terms 
of the two arbitrary constants cq and do, which are given 
by 

i 

J2 ai-fe[cfe(i - k) + dk] 
^ _ _k=i 

- 1 + Co) + do 
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J2 ai-k[ck{i - k - m) + dk] 

k=l 

{i — m){i — m — 1 + Co) + do ' 

dj(co — 1 — 2m + 2i) 
(j — m)(z — m — 1 + Co) + do 

i 

J2 {ai-feCfc + b,^k[ck{i - k-m,) + dk]} 

k=l 



{i ~ m){i — m — 1 + Co) + do 



(3.57) 



IV. FAILURE IN SOLAR SYSTEM TESTS 

The solar system tests are usually written in terms 
of the Eddington parameters, by following the so-called 
"parameterized post-Newtonian" (PPN) approach, intro- 
duced initially by Eddington [s^l- These parameters are 
often written in terms of the line element in its diagonal 
form. 



(4.1) 



Then, the gravitational field, produced by a point-like 
and motion-less particle with mass M, is given by 



„2* 



„2<I> 



1 + 27 



(4.2) 



where j3 and 7 are the Eddington parameters. For the 
solar system, we have Vg = GMq/c^ ~ 1.5 km, and its 
radius is tq ~ 1.392 x 10® km. So, within the solar sys- 
tem the dimensionless quantity x[= GM/ (rc^)] in most 
cases is much less than one, x ^ ''"g/lQ — 10^®. The 
Shapiro delay of the Cassini probe [55|], and the solar 
system ephemerides [11] yield, respectively, the bounds 

13, 



7 - 1 = (2.1 ±2.3) X 10"^ 
/?-l = (-4.1 ±7.8) X 10"^ 



(4.3) 



GR predicts = 1 = 7 precisely. To study the solar 
system tests in the HL theory, we may first transform 
the above experimental results in terms of the ADM line 
element with the projectability condition. 



ds^ = -df + e^'' {dr + e^-'^dty 



for which it can be shown that 144 



r'^dn'^. 



(4.4) 



1 



In <^ 2c2 



\ c^r 



(4.5) 



In the case A = 1, two different identifications were 
prescribed. One was to consider A as part of the metric 
via the replacement [s^ , 



(4.6) 



N ^ N -^A. 



With such an identification, the diagonal solution [381 144 



iV = 1, iV' = 0, / = 1 



2m 

r 



2to 
r 



A^\-Ao\ \ , ^ = 0, (A = l), (4.7) 



produces exactly the Schwarzschild solution in the form 



with * 



-$ = iln(/). Note that in writing 



Ea. (l4.7p . the speed of light appearing in Ea. (j4.6p had 
been set to one. As a result, the theory is consistent 
with observations [ssj . 

However, the solution (|4.7p is not unique, and there 
exists a larger class of non-diagonal solutions given by 
0, 



r = 



In /i(r) 

1 (2B 1 
2IH-+3 

0, = 0, (A 



Ar^ - 2A(r) 



A(r')dr' 



(4.8) 



where the gauge field A(r) is undetermined. If one does 
not consider the gauge field A as a part of metric [3] , but 
simply considers it as representing a degree of freedom of 
the gravitational field, as the Brans-Dicke scalar field in 
the Brans-Dicke theory [58], one finds that the above 
solutions are consistent with all the solar system tests, 
provided that |44| . 



A{r) = O 



(4.9) 



In the rest of this section we shall follow the second pre- 
scription, i.e., setting directly. 



(r, 17) = (//,!/), 



(4.10) 



and verify whether this prescription can be generalized 
to the case A ^ 1. As to be shown below, the answer is 
unfortunately negative. 

To this goal, we first note that the cosmological con- 
stant A has negligible effects within the solar system. In 
addition, the spatial curvature of the solar system is neg- 
ligible. In fact, for the metric (|4.4p . it takes the form, 



R = — 



1 



-2n 



(1 - 2rn') 



8(7-1)^ /gmv 

\ c^r / 



(4.11) 



for r >> To 



GM/c^. Note that in writing the last 



step of the above equation, we had used Ea. (j4.5|) . Thus, 
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in the solar system we have = R/2 < 10 km 
Therefore, without loss of generality, we set 

A = Ag = .g, ==0, (s>2). (4.12) 



Then, we obtain 



Then, the solutions are those given by Eqs. p.l6l) and 
p.34p . from which we find that 



1 / 2B\ 



+0 



(4.13) 



where e = Bc^ / [GM). Comparing the above with 
Ea. (li3)) . we find that 



7-1 



Be" 
GM 



< 10" 



for M — Mq. As a result, we have 

_ r _ 1 (GMy^ 
^^25" 2(7-1) \ ) ^ 

From the relations, 

F{a,b;d;x) = {I - x)-^F \b,d ~ a;d; 



F{a,b;d;l) = 



(4.14) 



(4.15) 



X- 1 



T{d)r{d-a-b) 

r(d- a)r(d- 6)' 



(4.16) 



where the last expression holds only for d ^ 
0, -1, -2, Re(d - a - 6) > 0, we find from Eq.dHJl) 
that h(r) has the asymptotical form. 



with 



1 

25 



h{r) ~ Dir, [x > 1), 
ai(-l)''i^(6,c- a;c;l) 



(4.17) 



+ a2(-l)''i^(6,c-a;c;l) 



7- Ao 



c=^.(4.18) 



Note that Ea. (|4.17l) can be also obtained directly from 
Eq.dS^Ml), which reads 

x^h![ + (a + 6 + \)xh!^ + afe/ii = 0, (4.19) 

for X ^ 1. Eq. (|4.19p has the general solution, 

hx{x) = d^x'^^'^^^l^ + rf2a;~^^+^°^/^ (4.20) 

where d\ and d2 are two integration constants. On the 
other hand, from Eq. p.29p we find that 



(4.21) 



h(x) = hQ{x)hi{x) ~ dix + 



dix, 



(4.22) 



which is precisely the solution given by Eq. (j4.17p with 
Di = di/{2B). Hence, we obtain 



r(r) 




(4.23) 



This is quite different from that given by Ea. (|4.5p with 
any choice of ai, a2, B, as long as A is not exactly equal 
to one. Therefore, the static vacuum solutions given by 
Eqs. (|3?T6)) and ((3T34)) with the condition (I4TT2)) is incon- 
sistent with the solar system tests, when the prescription 
is used. 



V. MOST GENERAL ANSATZ WITH 
SPHERICAL SYMMETRY 

In the previous section, based on the stationary ansatz 
(|3.1|) . we have seen that the prescription (|4.10l) leads to 
failure in solar system tests for A ^ 1, however small 
|A — 1| is. Hence, in the next section we shall consider 
another prescription. Before that, however, in this sec- 
tion let us consider the most general ansatz with spheri- 
cal symmetry and show that the prescription (|4.10l) never 
recovers the Schwarzschild geometry in the A — 1 limit 
with Ag = 0. This confirms that a prescription beyond 
(|4.10p is absolutely necessary. 

In order to find the most general ansatz, note that one 
can always choose time and spatial coordinates so that 
N — 1 and A^* = at least locally. One can also set 
(/? = by the U{1) gauge freedom. With 



N=l, Lp = 0, 



(5.1) 



it is obvious that the most general ansatz with spherical 
symmetry is 

g.jdx'dx^ = e2^(*'^'da;2 + e^'^^^'^'df^^, A = A{t,x). 

(5.2) 

Independent equations are the equation of motion for the 
gauge field A, the x-component of the momentum con- 
straint and the xx-componcnt of the dynamical equation. 



^ One must not confuse with the function B{t, x) used in this 
section and the constant B used in the expression of /(r) in the 
previous and next sections. 
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The equation of motion for the gauge field A is written 



as 



-2B+3C 



leading to the general solution 



F{t)+eC -{Ag/3)e^c 



0, (5.3) 



(5.4) 



where F{t) is an arbitrary function of time. The momen- 
tum constraint is 



d^dtC + dxCdtC ~ dxCdtB 



+ (A-1) 



0. (5.5) 



By using the solution (|5.4p . this equation is reduced to 
an equation for C: 



= (A - 1) ^ ci 
d.dtC 



+C2- 



dxC 



(dxCf [dxCf 



(5.6) 



where 



ci = ^ 



3F^e-^ + 2(3 - Age^^)F 



+ V^(A,e3C-3e^)2 



C2 



C3 = 



1 

6A 



1 

2A 
3 



63^26^*^ + 3(39 - llAge^C)^ 
3(3Age2'^ - l)F + 4Age3 



C4 - ^{l-Age^'')dtF, 
A = F + - (Ag/3)e 



3C 



(5.7) 



Finally, the dynamical equation can be considered as an 
equation determining the gauge field A. With the pre- 
scription (|4.10l) where A does not participate in the ge- 
ometry nor in solar system tests, the dynamical equation 
is not of our interest. 

Now let us expand C by (A — 1) as 



C(i,x) = ^C„(i,a;)(A-l)'^ 



(5.8) 



We shall see below that the momentum constraint equa- 
tion (|5.6p can be solved iteratively order by order in the 
(A — 1) expansion, under the condition (|5.12l) below. 



First, the zeroth order solution Cq is obtained as a 
solution to the following algebraic equation 



(5.9) 



where G{x) is an arbitrary function of x. Next, let us 
show by induction that the rt-th order solution C„(i,x) 
can be obtained by solving (|5.6p order by order. For this 
purpose, let us expand the expression inside the squared 
bracket on the left hand side of (|5.6I) by (A — 1) as 

oo 

Age^^ - 3F(t)e-^ = ^ e„(i,x)(A - 1)", (5.10) 

according to the expansion (j5.8p . It is easy to understand 
that Qn has the form 



Qn = {2Age^^'> + iF{t)e'^«) C„ -I- g„ 



(5.11) 



where Qn depends only on Ci (i = 1, 2, • • • , n — 1). Thus, 
provided that 

2Age'^^° -f ■iF{t)e-'^° ^ 0, 
F -I- e^^o - (Ag/3)e^<^'' ^ 0, 
dxC^ + 0, (5.12) 



we obtain 

Cn (ij 2;) 

Ai 



Ai 



Qn\C\, - ■ ■ , C„_i; i, x) 
d£ Sn\C\, ■ ■ ■ ,Cn^i;t' ,x) 



2Age^^° +3F{t)e-^", 



(5.13) 



where Sn is the n-th order part of the right hand side of 
(15. 6p . which also depends only on (i = 1, 2, • • • , n — 1), 
and to is an initial time. Note that the initial value of 
Cn at t = to has been set to zero by redefinition of G{x) 
and that the change due to the shift of the initial time 
to corresponds to redefinition of G{x). From this result, 
it is obvious by induction that the solution of the form 
(|5.8p can be obtained up to any order of the expansion. 

Let us consider the zero-th order solution (|5.9I) with 
Ag = 0. In order for the expansion w.r.t. (A — 1) to 
make sense, the condition (I5.12[) must be satisfied. In 
particular, F{t) should be non- vanishing. Otherwise, the 
denominator on the r.h.s. of (j5.13p would vanish. We 
thus assume that F{t) ^ 0. 

We would like to see if the Schwarzschild geometry is 
recovered in the limit A — ?> 1 with Ag = or not. One 
of the simplest ways is to calculate the 4-dimensional 
Einstein tensor for the 4-dimensional metric 

where 



27 F{tf{dxG{x)f 



G(x)4(3-G(x)) 



, Co = In 



G{x) 



(5.15) 
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Non-vanishing components of the 4-diniensional Einstein 
tensor are 



F2 



G 



(4)9 



1 

1 

F2 



{dtFf + 2FdlF 
{dtFf + 2FdlF 



27 

54 



(5.16) 



In order to recover the Schwarzschild metric, ah of these 
components must vanish, leading to dtF = G = 0. How- 
ever, in this case the regularity of Cq imphes that F = Q, 
contradicting with the assumption F ^ Q. Note that 
F ^ is a necessary condition for the continuity of the 
A — > 1 hmit. 

If we set _F = then the A ^ 1 limit is singular. Thus, 
the zero-th order solution is not obtained as the A — > 1 
limit of a solution with A 7^ 1. Instead, it represents 
a solution with exactly A = 1. In this case, (|5.9p with 
Kg — implies that G = and leaves Go unspecified. 
There is a choice of Go giving rise to the Schwarzschild 
metric. 

Just for completeness, let us consider the case with 
F = and 7^ 0. In this case the A — 1 limit is 
continuous. However, the zeroth order solution is 



A. General Coupling of A and tp to Metric 

To the above claim, we first note that under the Ufl) 
transformations, the ADM quantities transform as 38l . 

EH, 

5aN = 0, 5c,N^ = NVia, Sc^g^j = 0, 

5aA — a ~ N''\'ia, da^ = —a, (6-1) 

where SaF — F — F, a[— a{t, x)] is the generator of the 
local U{1) gauge symmetry. From the above we find that 



5aA = a-N'S/ia, 



where A is defined in Ea. (|1.6p . and 

a = A-A, JV' = N'+NW'ip. 



(6.2) 



(6.3) 



If we require that the line element ds^ be gauge- 
invariant not only under Diff(M, J^) (|1.2p . but also under 
the enlarged symmetry (|1.4p . then ds^ defined by. 



ds^ = ^APc^dt^ + g^j (dx' +Af'dt) {dx^ +Af^dt) 



has the desired properties, where 



(6.4) 



Bo 



1 



■In- 



3(a,G)2 



2 4AgG(G - 3) 



Co 



2 A„ 



resulting in 



G^f = -A 



9' 



(5.17) 



(5.18) 



Hence, the Schwarzschild metric is not recovered in this 
case unless the limit Ag — >■ is taken. If we take this limit 
then the A — >■ 1 limit becomes singular. Thus, again, the 
Schwarzschild solution is not obtained as the A — >■ 1 limit 
of a solution with A 7^ 1. Instead, it represents a solution 
with exactly A = 1. 

In summary, if A = 1 and Ag = exactly then the 
Schwarzschild metric is one of solutions. However, if we 
consider A 7^ 1 and Ag = 0, then the Schwarzschild metric 
is never recovered in the limit A — > 1. This conclusion is 
based on the prescription (|4.10p and thus implies that a 
prescription beyond (|4.10p is absolutely necessary. 



VI. PHYSICAL INTERPRETATION OF A AND 

In this section, we shall show that a proper general- 
ization of the prescription of (|4.6p can lead to solutions 
that are consistent with solar system tests even for A 7^ 1 . 
From such a generalization, the physical and geometri- 
cal interpretations of the gauge field A and Newtonian 
prepotential ip also become clear. 



J^=N~-{A-A)., 



(6.5) 



where v is a dimensionless coupling constant subjected 
possibly to radiative corrections. Similar to N, such de- 
fined Af is also dimensionless, [J\f] = 0. With this pre- 
scription, one can see that the Newtonian prepotential ip 
is tightly related to the shift vector A/"*, while the geo- 
metrical lapse function Af is related to both A and (p. In 
addition, since 

[dx] = -1, [dt] = -z, [c] = = z-1, 

[dt] 

[N] - 0, [TV'] =z-i, [.g.j] = 0, 

[A] = [A] = 2{z-l), M=z-2, 
\a\ = z-2. 



we find that 



[ds] = -1, 



(6.6) 



(6.7) 



i.e., it has the dimension of length. Moreover, with the 
gauge choice 1^ = and setting v — 1, Eq. (|6.5p reduces 
to Eq.gSl). 

In the Newtonian limit, we have [3a ^31 



500 



(6.8) 



in the coordinates x^ = {ct, a;*), where e = \v/c\ <C 1, and 
V denotes the typical velocity of the system concerned. 
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Comparing it with the metric given by Eas. (|6.3p - (|6.5I) . 
we find that the Newtonian potential cj) is given by 



with 



-v{A + ^)~ -N'N, 
+ {v-l) (n' + ivV ) V,(p, (6.9) 



A^ = l, -\2N, + W.M = 0{e). (6.10) 

c 

To study further the meaning of the above prescription 
and the physical interpretations of the gauge field A and 
Newtonian prepotential f, let us turn to the solar system 
tests again. 



B. Solutions with the Gauge ^ = 

For the ADM decomposition p.ip without fixing the 
U(l) gauge, there are three independent equations, given 
by Eas. (|D.l| ) - (|D.3p . To solve these equations, let us 
first note that the prescriptions of Eqs. (l6.3p - (|6.5p do 
not change the spatial metric gij. As a result, the con- 
straint on the spatial curvature R takes the same form 
of Eq. (14.111) . Therefore, in the present case the condition 
(|4.12p can be still imposed safely. In particular, with the 
gauge A — [cf. Appendix C for different gauge choices.], 
Eqs. (|D.ip and (|D.2p for A = 1 have the solutions. 



2B 



/-I , h^-fip', (A = l, A^O), 



.11) 



where ip must satisfy the dynamical equation Egs. ljD.Sp . 
which now reads. 



3.12) 



The general solutions are given by, 
2r 



ip{r) ^ipo± 



r-2B 



r-2B 



1/2 



dr, (6.13) 



where ipo and (pi are integrations constants. Substituting 
the above into Ea. (|0| we find that TV* = 0. Then, 
Eq. ()6.4p reduces to. 



dr^ 



r^d^Vt, 



where 



I 



fir 

, (A = l, A = 0). 



3.14) 



2- w I 2 
2B 



3.15) 



When X ^ 1, we have 



8T2 



J = l + ex + eV + 0(x') 



(6.16) 



where T = v{l + fi/2) — 1. The factor T appearing in 
the expression of Af can be dropped by rescaling t — Tt. 
Then, comparing Ea. (|6.16l) with Ea. (|4.2p we find that 



B = -f 



fGM\ 
2 



2+(7-l)v?i' 



3.17) 



For ifi ~ we obtain the constraint |u — 1| < 

O (lO^^) from (|4.3p . For extremely large value of tpi, 
say <pi ~ C (10^), \v-l\~0 (1) is also allowed. How- 
ever, we consider this large value of fi unrealistic and 
consider the case with ipi — O (1) only. Note that the 
Schwarzschild solution corresponds to B = GM /c? , v — 



C. Solutions with the Gauge tp ■ Q 

On the other hand, in the case A ^ 1 let us consider 
the gauge /i = 0. Then, Eqs.(|D.ip and (|D.2p yield. 



f{r) = 1- 



2B 



ipir) — I r ^ 



-bo r ^ F 



, „ , 9 — Ao — 3 — Aq 2 — Aq 
Ol F I : — ; : : : ; X 



dr 



4 ' 4 2 
M^/3 + Ao 9 + Ao 2 + Ao 



(6.18) 



where 61^2 are constants, and Ao is given by Eg. (13. 331) . 
Substituting the above into Ea. (|D.3p . we find that 



^(0 = ^/1-^ Mo - 



P(r) 



2B 



--dr , (6.19) 



with 
P(r) 



(21 -33A)B^ + 2(3-4A)r 



4(2B- r)r2 
-^2Sr(16A - 11)1 {if'f ~ 4:Br 

-r^ir-2Bfi\~lwA 



-t-^^-^^' \(2B - r)rp"' - 2Bp"] 
2r 



(6.20) 
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When bi 



0, the above solutions reduce to 



ip{r)^ipo, A=l-Ao\/l-^. 



3.21) 



Substituting it into the metric (j6.4l) , and considering the 
gauge choice ft, = 0, we find that it takes exactly the 
form of Eq. (|6.14p with the metric coefficients given by 
Eq. ()6.15|) and with the replacement (pi — > —2Ao. Thus, 
the PPN parameters (3 and 7 are given by (|6.17p with (pi 
replaced by —2Ao. For Aq ~ C(l), we obtain the con- 
straint |u - 1| < O (10~^) again from (|4.3|) . Therefore, 
the prescription (|6.4|) leads to consistent results with so- 
lar system tests even for A 7^ 1. 



VII. CONCLUSIONS 

In this paper, we have studied spherically symmetric, 
stationary vacuum configurations in the general covari- 
ant theory of the Hoi^ava-Lifshitz gravity with the pro- 
jectability condition N = N(t), and an arbitrary value 
of the coupling constant A |38l-l4l|. In particular, in Sec. 
Ill we have obtained all the solutions with the assumed 
symmetry in closed forms. 

When applying these solutions to the solar system tests 
(Sec. IV), we have shown explicitly that the ADM- type 
identification (|4.10|) between the metric coefficients and 
the basic quantities N, and do not render the A 7^ 1 
solutions consistent with solar system tests, no matter 
how small |A — 1| is. (On the other hand, when A = 1 
exactly, there is a spherically-symmetric, stationary vac- 
uum solution which is consistent with the solar system 
tests [3.) 

To show that this is indeed the case in more general 
situations, we have devoted Sec. V to consider the most 
general ansatz (|5.ip and (|5.2p with spherical symmetry 
and shown that the prescription (14.101) never recovers the 
Schwarzschild geometry in the A — > 1 limit with Ag = 0. 
Thus, one needs either to invent a mechanism to restrict 
A precisely to its relativistic value Xgr = 1, or to consider 
the gauge field A and/or the Newtonian prepotential ip as 
parts of the 4-dimensional metric on which matter fields 
propagate. 

In the case A = 1, HMT proposed the identification 
(|4.6p jssj but clearly it is not gauge-invariant. Requir- 
ing the line element be gauge-invariant not only under 
DifF(Af, J^) (|1.2p . but also under the enlarged symmetry 
(II. 4p . in Sec. VI we have proposed the identification 
(16. 4p . where w is a dimensionless constant to be con- 
strained by observations/experiments. When w = 1, it 
reduces to (14. 6p in the gauge ip = 0. Applying such a pre- 
scription to the cases both with A = 1 and with A 7^ 1 , we 
have shown that the resulted metric is indeed consistent 
with the solar system tests, provided that |w— 1| < 10~^. 
With such identifications, one can also see the physical 
and geometrical roles that A and ip play. In particular, 
the Newtonian prepotential p is tightly related to the 



shift vector, while the geometrical lapse function is 
related to both A and p. 

Finally, we note that it still remains to be understood 
how to obtain the prescription (|6.4p (with w ~ 1) from 
the action principle ^ . Actually, in the UV, N and A — A 
have different scaling dimensions and thus, it is not easy 
to imagine how their linear combination can universally 
enter the UV action of matter fields. On the other hand, 
in the IR, N and A — A have the same scaling dimensions 
(they are actually dimensionless) and thus, the prescrip- 
tion (j6.4l) is not forbidden a priori. It is therefore im- 
portant to investigate whether the prescription (j6.4p can 
emerge in the IR and, if it does, how. 
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Appendix A: Field Equations 

Corresponding to the actions (|2.5p . the Hamiltonian 
and momentum constraints are given respectively by. 



d^x^g [Ck + Cv~ pG'^V^VjP - (1 - a) (vV)^ 
= 87rG j d^xy/ljj\ (A.l) 



where 



T^i] - ^Qi] - (1 - V^(/3 



(A.2) 



J* 



J,. = -N 



6N ' 
Kij + XKgij, 

SCm 



(A.3) 



Variation of the action 
yield, respectively. 



with respect to p and A 



g'^ + v,v,^) + (1 - A) V2 (^K + 



= SttG J. 



i?- 2A„ = SttGJa, 



(A.4) 
(A.5) 



In [6011 the coupling of the HL covariant theory with matter was 
considered from the action principle. It was shown that Newto- 
nian gravity cannot be recovered in the weak gravitational field 
approximation, based on several assumptions, including the one 
that the coupling among matter, the gauge field A and the New- 
tonian prepotential tp be described by the recipe provided in [40ll . 
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where 



J.^-^. J. ^2^1^. (A.6) 



Sip ' "^-^ 5A 
On the other hand, the dynamical equations now read 



Appendix B: G and H defined in Eas. (|3.14[) and 



The functions G and H defined in Eas. ljXTI)) and 
are given by 



V5 TT^J' - vQ'^ - (1 - A)g*^ VV 



= -2 (i^^yj + 2\KK'^ 



N 



-2(1 -A) (/f + VV)V*VV + ^'^V2(^ 



+ (1-A) 

+ F'' + f;J' + F^^' + 8^Gr*^' , (A.7) 



where [K^Y' = K'^Kj, = + /2, and 



V9 -^ffy t^n 



(ip,n)' 



n=l 



i<r + vV)vv + — v>, 



1 

N 



(A.8) 



with = (2,0,-2,-2,-4,-4,-4,-4,-4). The stress 
3-tensor r*-' is defined as 



V9 



(A.9) 



and the geometric 3-tensors (F^) and F^*^ are given 
in m. 

The matter components (J*, J', Ji^, J^, r*^) satisfy 
the conservation laws, 



= 0, (A.IO) 



-^Vfej'= + J^V,(/.-^V,A = 0. (A.ll) 



G(r) 



r/i[2/i' - r(A - l)/i"] + -r^{X - l)/i'' 



r 

2/ 


r(A 




8/2 ^ 


1 




4(4 


r4 j 



8(46.g4 + 1755 + 7ge + 28<?7 + Sgs)/^ 



-4[8gy/Wr4 + S^s/^r^ + Ug,C\^ + Bg^C'r^ 

+(48ff4 + 14^5 + 3^6 - 4857 - l(igs)f"r^ 
-2(4854 + 1455 + 356 - 4857 - 1658)/V + I8O54 
+6655 + 2756 + 4857 + 1258] /' + [ - 3(857 
+358) (/")''^' + (4854 + 2255 + 1256 + 4857 
+1358)(/')''^' + 4(852r2c2 + 353r2c2 + 4854 
+1455 + 356 - I657 - 458)/"r2 - 2/' [3(3254 
+1255 + 556 - 58)/'V2 ^ 2((857 + 358)/"V3 
+9654 + 2855 + 656 - 3257 - 858)]r + 4[r4c4 
+1252r2c2 + 453r2c2 + 8454 + 3O55 + I256 - 857 
-658]]/ + (3254 + 1255 + 556 - 58)r'(/')' 
+4[C^Ar6 - + 252(^2 + 53CV2 + 454 
+255 + .96] + r'if'f [ - 852C''^' - 353C' 
+(857 + 358)/'V2 - 4854 - 1455 



-356 + 1657 + 453] I , 
H{r) = 2Ar2 /i/i" + (A + l)r2 (/i')^ +4 (2A- l)r/i/i' 



(2A + l) r2f 
/ 



hh' 



\rf" + (2A - 1)/' 



+ 



32(4654 + 1755 + 756 + 2857 + §58)/' 



4^4^4 

+4[857/(5V5 + 358/(5V5 + 4854/'%^ 



356/"' - 



+1455/ 
-1658/"' 

+1455 + 356 - 4857 - 1658)/" 



'/ r-3 



2C'r2 



- 4(4854 
6(254 - 355 
-456 - 7657 - 2558)/V - 72O54 - 26455 
-10856 - 19257 - 4858] /' + 2 [ - 1652C'/'"r^ 
-653C'/'"r' + 3(3254 + 1255 + 556 - g&){f"fr^ 
-Q&gifr^ - 2855/'" - 656/'%^ -I- 3257/"' 
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+8ff8/"V' + 48g2C'^' + 1653C'r-2 - 3(112g4 
+3035 + 4.96 - 144.97 - 47gs)ifTr^ + /"(5(8g7 
+3.98)/'" + 8(4854 + 14.95 + ^96 ~ 16g7 

-^98)y + f ^3(32.94 + 12.95 + 5.96 - gs)f"'r^ 

+(4834 - 2.95 - 1536 - 240^7 - 7438)/"^' 
+2(3(837 + 338)/(4V4 + 2(14.92r2c2 + 553^2^2 

+36.94 + 24.95 + 18.96 + 9637 + 2438))^ r + 67234 

+240.95 + 96.96 - 6437 - 4838] / - (I6.94 + IO35 
+736 + 48.97 + 14.98)r3(/')3 + 8{~C^Ar^ + 232C'r 
+93C^r^ + 8.94 + 435 + 2.96) - r/' [ - (837 
+338)(/")''^' + 2(8.92r2c2 + Sg^r^C^ + 4834 
+1435 + 336 - I637 - 438)/"r2 + iir'^C" 
+12g2r\^ + 433^2(2 + 8434 + 3O35 + I236 

„3 



When h = 0, we have 



a{t,r) 



Ait',r)dt' + aoir), (C.4) 



where ao{r) is an arbitrary function of its indicated ar- 
gument. Thus, in this case the gauge is fixed only up to 
an arbitrary function of r. 

When h ^ 0, we introduce two new variables u and v 
via the relations, 



dt = Gdv + Fdu, 
dr = h{Gdv-Fdu), 



(C.5) 



where F and G are functions of u and v only, and satisfy 
the integrability conditions, 



F^v — G,u = 0, 
(F/i),„ + (G/i),„ = 0. 



(C.6) 
(C.7) 



837 - 638)] + 3r2(/')2 [(837 + 338)/"V' 
+(3234 + 1235 + 536 - 38)/'V' + 9634 + 2835 

+636 - 3237 - 838] I, (B.l) 



Note that one should not consider Ea. (|C.5[) as coor- 
dinate transformations, because they are forbidden by 
DiS{M,F), but rather a technique to solve Eq. (jC.3p . 
Then, in terms of u and v, Ea. (IC.3p takes the form, 
a u = —FA, which has the solution. 



where /(") = 



Appendix C: The U(l) Gauge Transformations and 
Gauge Choices 

Under the U(l) gauge transformations (|6.ip . 
in the spherically symmetric case, the variables 
{N, iV*, gij , A, (f) transform as, 

Sa^N = 0, Sa^N' = Slfa', 5^g,j = 0, 

SaA — a — ha', Sa'^ = —a, (C.l) 

where a = a{t,r). From these expressions, one can see 
that various gauges can be chosen. 



A. = 

In this gauge, we have 

a = v{t,r), (C.2) 
which is unique, and is the gauge used in Section III. 



a{t,r) 



F{u',v)A{u',v)du' + ai{v), (C 



where ai is an arbitrary function of v only, and u — 
u{t,r) and v ~ v{t,r), given through Eqs. (IC.5P - (IC.7p . 
Therefore, in the present case the gauge is fixed up to an 
arbitrary function of v. 



C. h^O 



In this gauge, we have 

h 

which has the solution, 

r h{ty)dr' 



a{t,r) 



f{ty) 



(C.9) 



(C.IO) 



where 02(1) is an arbitrary function of t only. 



Appendix D: Field Equations without Specifying the 
U(l) Gauge 



A — u It can be shown that in the spherically symmetric case, 

there are only three independent field equations: the con- 
In this gauge, we have straint obtained from the variation of the gauge field A 

given by Ea. (IA.4[) . the momentum constraint (jA.2p . and 
a — ha' — —A. (C.3) the rr-componet of the dynamical equations (|A.7p . For 
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the ADM decomposition given by Ea. p.l[) . they read, 
respectively, 



(r/)'-(l-A,r2) =0, 



(D.l) 



il-X)lr^fh"-'-^{rf'-4f)h'- 



2f 



^y(/')^ + y/.f 



Af^' - 4r/' 
r-V"-(4r/ + 3rV')^"-2r/(p"' 



-rff'h + f if + = 0, 

r r^j 



(D.2) 



o!,2 



rf 

(4654 + 1755 + 756 



'J7 + 9gs)f - 4( - (857 - 358)/(*V4 + (1452 

+553)C^?-^ + (4854 + 1455 + 356 + 4857 
-1658)(r/" - 2/')r + I8O54 + 6655 + 2756 

+ 12(58 - 457)) f + {i{r*C' + 4(352 + 93^^ 

+8454 + 3O55 + 1256 + 857 - 653) + r((4854 



+2255 + 1256 - 4857 + 1358)r(/')' - 2(2(358 
-857)/"V3 + 3(3254 + 1255 + 556 - g8)f"r^ 
+4(4854 + 1455 + 356 + I657 - 458))/' 
+r/" (3(857 - 358)/"r2 + 4((852 + SgsV^C^ + 4854 
+1455 + 356 + I657 - 458))))/ + 4(r2(r2(r2A - 

+252 + 33) C + 454 + 255 + 36) - iff ((852 
+353)r-^C^ + 4854 + 1455 + 356 + I657 - 458 



+r((-3254 - 1255 - 556 + 58)/' + (857 - 358)r/")) 
„ iff , 8/' 32/ 



+ (1-A) 



4/" 



/ 



r 



+ 8/VV" 



^2J3 



3r2(/')2-32/2-4r/(r/"-2/') 



^hh"^f^^f' 64 6(/')^\ 



-8(/i + /^')^" 



0. 



(D.3) 
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